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ALS deep UR. 
(云南 农业 大 学 水 利 学 院 ,650201 昆明 ) 


摘 要 :为 了 研究 三 维 串联 多 圆柱 绕 流 在 不 同 圆柱 间距 比 L/D(L 为 圆柱 的 间距 ,D 为 圆柱 直径 ) 情 
况 下 的 尾 涡 结构 ,基于 格子 玻 尔 兹 曼 (LBM) 方 法 ,采用 大 涡 模 拟 (LES) 方 法 计算 了 高 雷诺 数 下 三 维 
单 圆柱 及 三 维 串 联 多 圆柱 在 不 同 L/D 情况 下 的 尾 涡 流 场 。 数 值 结果 表明 :LBM-LES 方法 可 以 较 好 
地 模拟 出 三 维 单 圆柱 以 及 串联 多 圆柱 的 尾 涡 结构 ;高 雷诺 数 下 单 圆柱 绕 流 尾 迹 变化 具有 不 稳定 性 ; 
串联 多 圆柱 绕 流 存在 临界 间距 比 L/D =4, 当 圆柱 间距 比 小 于 临界 间距 比 时 ,圆柱 间 不 会 产生 尾 涡 
脱落 ; 当 圆 柱 间 距 比 大 于 临界 间距 比 时 才 会 产生 明显 的 涡 脱 落 现 象 , 且 前 方圆 柱 脱 落 的 涡 对 后 方圆 
柱 有 冲击 作用 ,会 对 后 方圆 柱 的 升 / 阻 力 系数 振动 幅 值 有 一 定 的 影响 ;串联 多 圆柱 绕 流 的 上 下 游 时 
均 阻 力 系数 曲线 在 临界 间距 比 时 会 发 生 转折 , 呈 和 逐渐 增 大 趋势 ; 当 圆 柱 间 距 比 小 于 临界 间距 比 时 ， 
在 下 游 圆 柱 前 增加 圆柱 会 对 下 游 圆柱 以 及 中 游 圆柱 的 时 均 阻 力 系数 产生 影响 ,而 当 圆 柱 间距 比 大 
于 临界 间距 比 时 ,在 下 游 圆 柱 前 增加 圆柱 对 下 游 圆柱 以 及 中 游 圆柱 的 时 均 阻 力 系数 基本 没有 影响 。 
关键 词 :格子 玻 尔 效 曼 方法 ;大 涡 模 拟 ;圆柱 绕 流 ; 尾 涡 结构 ;数值 模拟 
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Research on wake structure of three-dimensional tandem 
multi-cylinder flow based on LBM-LES 
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Abstract :In order to study the wake structure of three-dimensional tandem multi-cylinder flow at different 
L/D (L is the spacing between the cylinders and D is the diameter of the cylinder) , based on the lattice 
Boltzmann method ( LBM) ,the large eddy simulation ( LES) method was used to calculate the wake flow 
field of three-dimensional single cylinder and three-dimensional series multi-cylinder under different L/D 
conditions. The results show that LBM-LES method can well simulate the wake structure of three-dimen- 
sional single cylinder and multi-cylinder in series. The wake change of single cylinder with high Reynolds 
number is unstable. There is a critical spacing ratio L/D 24. When the spacing ratio between cylinders is 
smaller than the critical spacing ratio, there is no vortex shedding between the cylinders. When the critical 
spacing ratio L/D is 4 ,the vortex shedding phenomenon is obvious , and the vortex falling off the front cyl- 


inder affects the rear cylinder, which will have a certain impact on the vibration amplitude of the lift and 
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drag coefficients of the rear cylinder. When the ratio of cylinder spacing is smaller than the critical spac- 


ing ratio ,the time average resistance coefficient curve of the upstream and the downstream of the flow a- 


round multiple cylinders in series will turn at the critical spacing ratio. When the boundary spacing ratio is 


4 , adding a cylinder in front of the downstream cylinder will affect the time average resistance coefficient of 


the downstream cylinder and the midstream cylinder. But when the cylinder spacing ratio is greater than 


the critical spacing ratio of 4 , adding a cylinder in front of the downstream cylinder does not affect the time 


average resistance coefficient of the downstream cylinder or the middle cylinder. The research results have 


certain theoretical reference value for practical engineering. 


Key words:lattice Boltzmann method; large eddy simulation; flow around a cylinder; wake vortex struc- 


ture ; numerical simulation 


多 柱 绕 流 在 实际 工程 中 有 较 多 的 应 用 ,例如 : 桥 
墩 、 海 上 石油 平台 支柱 及 一 些 桩 基 式 码头 。 水 流 在 
流 过 多 圆柱 结构 时 会 在 圆柱 后 产生 交替 脱落 的 波 
沽 .周期 性 脱落 的 流 涡 会 产生 脉动 压力 ,可 能 会 使 结 
构 物 发 生 大 幅度 的 振动 ,严重 情况 下 可 能 会 破坏 束 
体 芋 筑 。 另 外 ,多 柱 绕 流 特性 与 单 圆 柱 绕 流 具有 明 
丑态 同 。 因 此 ,研究 多 圆柱 绕 流 问题 对 实际 工程 具 
KREWE. 

CD 目前 ,国内 外 研究 人 员 大 多 通过 物理 模型 试验 
和 和 北 值 模拟 的 方法 来 研究 单 圆 柱 绕 流 和 多 圆柱 绕 
流 N 对 于 物理 模型 试验 的 研究 , 杨 群 等 "运用 刚性 
模型 测 压 风 洞 试验 方法 对 单 圆 柱 、 不 同 间距 下 串 列 
双 圆 柱 和 串 列 三 圆柱 绕 流 的 时 均 压 力 分 布 和 气动 力 
进 各 了 研究 。 研 究 发 现 不 同 间距 下 存在 一 个 临界 间 
距 痢 界 间距 前 后 的 流 态 完全 不 同 , 且 时 均 压力 分 布 
和 时 均 阻力 系数 有 着 明显 的 不 同 。 沈 国 辉 等 ”通过 
在 天 匀 流 中 进行 不 同 风速 下 的 测 压 风 洞 试验 来 研究 


光 清 圆柱 的 气动 力 系数 和 绕 流 特性 。 人 研究 发 现 层 流 


涡 脱 是 规则 的 ,而 汕 流 涡 脱 是 随机 的 。LOPES 等 中 
通过 固定 试验 与 数值 模拟 结合 对 多 柱 构 型 进行 研 
究 。 结 果 表 明 多 柱 构 型 没有 规则 的 涡 脱 落 出 现 , 且 
柱 截面 类 型 决定 涡 脱落 行为 。 除 了 运用 物理 试验 
外 ,许多 学 者 也 会 运用 数值 模拟 方法 来 研究 圆柱 绕 
流 问 题 ,并 且 也 都 得 到 了 一 些 有 参考 价值 的 成 果 。 
曹 兴 等 所 采用 数值 模拟 的 方法 研究 了 脉动 流下 圆柱 
绕 流水 动力 特性 。 结 果 表 明 改 变 脉 动 流 频率 可 以 对 
涡 脱 落 速度 有 影响 , 且 脉 动 频 率 跟 量 纲 为 一 的 脉动 
振幅 的 增加 会 使 升 /阻力 振幅 增 大 。 田 启 龙门 采用 
离散 涡 方法 对 二 维 双 圆 柱 水 动力 特性 进行 数值 模 
拟 。 结 果 表明 间隙 比 对 圆柱 的 时 均 阻 力 大 小 有 影 
响 , 间 际 比 过 大 时 圆柱 干涉 将 变 小 。Abbasi 4&5 35 
用 数值 模拟 的 方法 对 串 列 圆柱 进行 研究 发 现 串 列 贺 


柱 绕 流 存 在 一 个 临界 间距 值 会 改变 流体 的 流动 特 
性 。 大 量 的 研究 结果 表明 ,数值 模拟 结果 与 物理 试 
验 结果 比较 接近 , 且 费 用 代价 较 小 ,说 明 运用 数值 模 
拟 的 方法 研究 圆柱 绕 流 问题 是 可 行 的 。 

早期 的 圆柱 绕 流 研究 多 以 模型 试验 为 主 。 由 于 
科技 的 快速 发 展 , 数 值 模拟 成 为 研究 圆柱 绕 流 问 题 
的 主要 手段 。 以 往 大 多 数 研 究 均 仅 仅 运 用 二 维 模型 
研究 圆柱 绕 流 问题 ,但 三 维 圆柱 绕 流 在 尾 涡 结构 及 
形态 上 与 二 维 圆柱 绕 流 存在 差别 。 近 年 来 ,许多 学 
者 开始 运用 数值 模拟 方法 研究 三 维 圆柱 绕 流 问题 。 
3E x ^" 38 JH OpenFOAM 模拟 研究 不 同 雷 诺 数 
下 三 维 圆柱 绕 流 特性 ,发 现 圆柱 的 展 向 流速 变化 引 
起 的 三 维 效应 使 得 升 /阻力 系数 小 于 二 维 升 /阻力 系 
数 。 涂 佳 黄 等 中 运用 大 涡 模 拟 方法 研究 亚 临 界 雷 诺 
数 三 维 串 联 圆柱 绕 流 。 研 究 表明 间距 比 的 变化 会 使 
双 圆 柱 结构 表面 流体 力 系数 时 间 和 空间 的 变化 规律 
性 慢 慢 减弱 ,而 临界 间距 比 时 流体 力 系数 变化 会 有 
较 强 的 规律 性 。 除 此 之 外 ,文献 [9-10] 运 用 大 涡 模 
拟 方法 对 三 维 圆柱 绕 流 问 题 进行 研究 ,结果 表明 大 
涡 模 拟 更 能 体现 出 细致 的 流 场 结 构 。 大 量 研究 表明 
三 维 圆柱 绕 流 模拟 明显 比 二 维 更 符合 实际 。 

国内 外 学 者 在 研究 圆柱 绕 流 问题 上 多 采用 
Fluent 等 软件 进行 模拟 探究 ,采用 LBM-LES 方法 研 
究 圆 柱 绕 流 的 文献 还 较 少 。 格 子 玻 尔 兹 曼 (LBM) 方 
法 是 近 十 几 年 来 国际 上 发 展 的 一 种 流体 系统 建 模 和 
模拟 新 方法 ,与 SPH .MPS 法 同样 采用 无 网 格 计算 方 
法 。 这 大 大 提高 了 工作 效率 , 且 该 方法 是 介 于 流体 
的 微观 分 子 动力 学 模型 和 宏观 连续 模型 之 间 的 介 观 
模型 , 兼 具 二 者 的 优点 ,克服 了 早期 的 格子 气 自 动机 
的 一 些 不 足 , 具 有 许多 常规 方法 无 法 比拟 的 独特 优 
势 。 周 凯 等 "1 运用 玻 尔 兹 曼 方法 对 二 维 串 联 圆 柱 
绕 流 进行 数值 模拟 研究 ,得 到 其 升 /阻力 系数 尾 涡 及 
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流 线 变化 ,结果 表明 串联 圆柱 存在 临 
"gsx in FE A 2 88 7 ONE EE VU 
研究 。 研 究 表明 不 同 间距 比 下 四 圆柱 


ABB. deu 
圆柱 扰 流 进行 
的 绕 流 形态 也 


不 同 , 且 存 在 一 个 间距 比 使 得 各 圆柱 的 阻力 系数 最 
小 。Zhang 等 ”运用 玻 尔 效 曼 方 法 对 可 渗透 菱形 柱 
体 非 定常 流动 进行 数值 模拟 研究 。 结 采 表 明 增 加 渗 


1 数值 模拟 方法 


1.1 格子 玻 尔 兹 曼 方 法 
本 文 三 维 圆 柱 绕 流 模拟 采用 三 维 19 速度 不 可 
压缩 格子 Boltzmann 模型 (D3Q19 模型 ) ,如 图 1 所 


透 率 会 抑制 涡 脱落 , 且 玻 尔 兹 曼 方法 模拟 结果 与 已 
有 数据 对 比较 吻合 ,说 明 该 方法 用 来 研究 圆柱 绕 流 
问题 是 可 行 的 。 
本 文 将 在 格子 玻 尔 效 曼 方法 (LBM ) 基础 上 , 结 
合 大 涡 模 拟 (LES ) 对 三 维 单 圆 柱 和 三 维 串 联 多 圆柱 
在 不 同 L/D 情况 下 的 绕 流 进行 数值 模拟 ,得 到 其 
升 / 阻 力 特性 及 尾 涡 结构 特征 。 


ed 
1 
0 
到 1 D3Q19 模型 
: Fig.1 D3Q19 model 
COE bs 4 PROC 
sq e;u (e;u)* u 
f=pollrost cr z © 


AP: w, 为 权 系 数 , wo z 1⁄3, w6 m 1/18 , w719 - 
1/36; C, =c/V3 为 格子 的 声速 ;流体 宏观 密度 p 和 速 
HE u 表示 为 


示 。LBM 模型 的 演化 方程 为 049 

f, (x t e; Nt t AD) -f(x,t) 

= 二 LACcD 2680] (1) 
式 中 :为 粒子 速度 分 布 函数 这 ”为 局 部 平衡 函数 ;r 为 
量 纲 为 一 的 松弛 时 间 ;e 为 离散 速度 ;At 为 时 间 步 长 。 
三 维 模型 D3Q19 模型 包含 19 个 速度 矢量 , 即 
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式 中 :4r = CLA, 为 过 滤 尺 度 , C, 为 WALE 模型 常 
数 ,通常 为 0.325 ,A, 为 单位 网 格 尺度 ; Gs 为 速度 梯 
度 张 量 平方 的 痕迹 对 称 部 分 ;Sy 为 分 解 尺度 的 应 变 率 
IKE; pI Kronedcker ò 函数 ;wy 为 张 量 下 标 。 


2 数值 计算 模型 与 方法 验证 


2.1 数值 计算 模型 
本 文 模拟 选 定 的 计算 域 尺 寸 为 4mxlmx2m， 
双 圆 柱 绕 流 计算 域 如 图 2 所 示 ( 单 圆柱 绕 流 和 三 圆 


柱 绕 流 计算 域 与 之 类 似 )。 图 


中 , 为 双 圆 柱 的 间 


X. 
Y 6f" 2 pu 


(4) 
(5) 


1.2 imb ES 

Am Uit E A HIIS UL C LES) , 该 方法 引 

入 了 一 个 附加 的 黏度 , 称 为 清流 涡 黏 度 ,以 模拟 亚 

格子 应 流 。LES 方案 是 Wall-Adapting 壁面 自 适应 局 

部 涡流 (WALE ) 医 度 模型 ,提供 了 一 致 的 局 部 涡流 
黏度 和 近 壁 行为 ,具体 方程 如 下 。 

， (G! G! jos 
V, =A; GU. (GG. 


(6) 


距 , 圆 柱 的 直径 也 =0.08 m, 圆 柱 的 高 度 为 2 m, A 
口 的 来 流速 度 =6 mss WENEN , RE u = 
0.016 Pa . s ,密度 p =1 000 kg/m 。 雷 诺 数 Re =3 x10'。 
圆柱 中 心 离 和 人口 距离 为 1 m, 离 上 下 边界 各 0.5m, 
边界 条 件 如 下 :入口 边界 条 件 为 速度 入 口 ,出 口 边界 
条 件 为 压力 出 口 , 侧 向 边界 为 对 称 边界 条 件 ,圆柱 表 
面 是 固定 无 滑 移 的 壁面 条 件 。 为 了 改善 近 壁 面 区 的 
清流 渐 近 行为 ,采用 增强 壁面 郴 数 对 圆柱 近 壁 面 进 
行 修正 。 对 于 串联 多 圆柱 绕 流 模 拟 取 2 .2.5 3、3. 5、 
4,5.6,7,8.9,10,3t 11 © L/D 间距 比 来 进行 模拟 。 


ChinaXiv 合 作 期 刊 go 


t=0.14s 时 ,圆柱 后 开始 慢 慢 生成 大 尺度 的 涡 结构 ， 


图 2 双 圆 柱 绕 流 计 算 域 


Fig.2 Calculation domain of double cylinder flow around 
2.2 方法 验证 

> 采 用 与 文献 [1] 相同 条 件 模拟 单 圆柱 和 串联 双 
1] 进行 对 比 。 单 


圆 社 模拟 的 阻力 系数 为 1. 15, 文 献 [1] 的 结果 为 
1O ,说 明 在 单 圆柱 情况 下 模拟 结果 与 试验 较为 接 
过 串联 双 圆柱 采用 2.2.5.3.3.5.4.4.5.5.6.7.8、 
1042 , 共 12 个 不 同 的 LD 间距 模型 进行 模拟 验证 ， 
模 徇 结果 如 图 3 所 示 。 可 以 看 出 串联 双 圆柱 在 不 同 
I-A 下 得 到 其 时 均 阻 力 值 ,利用 LBM-LES 方法 模 
WERSU 1] 进行 对 比 发 现时 均 阻力 系数 曲线 
变 北欧 势 基本 接近 ,可 以 得 出 本 文采 用 的 方法 是 可 
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图 3 串联 双 圆 柱 模拟 结果 


Fig.3 Simulation results of tandem double cylinders 


3 数值 模拟 结果 及 分 析 


3.1 三 维 单 圆柱 绕 流 结果 
图 4~ 图 6 分别 为 单 圆柱 绕 流 典型 时 刻 流向 截 
面 和 展 向 截面 的 瞬时 涡 结构 分 布 图 。 可 以 看 出 ,在 


随 着 时 间 的 推移 , 尾 涡 逐 渐 向 后 发 展 ;在 1=0.22s 
时 ,圆柱 下 游 大 尺度 涡 结构 逐渐 破碎 形成 小 尺度 的 
涡 ; 在 1=1.24 s 时 , 单 圆柱 尾 涡 结构 已 经 向 下 游 区 
域 充 分 发 展 ,形成 类 似 卡 门 涡 街 , 尾 迹 结构 变 得 不 稳 
定 , 可 能 会 引起 涡 激 振动 。 


«Qv 


Ca) 流向 截面 


(b) 展 向 截面 
图 4 1:20.14 s IAIN 


Fig.4 Vortex structure distribution at t 20.14 s 
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Ca) 流向 截面 


(b) 展 向 截面 
5 120.22 s 时 涡 结 构 分 布 


Fig.5 Vortex structure distribution at t 20.22 s 
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(b) 展 向 截面 
图 6 1=1.24s 涡 结构 分 布 
Fig.6 Vortex structure distribution at t = 1.24 s 
3 三 维 串 联 多 圆柱 绕 流 结果 
@ 〇 通过 模拟 不 同 间距 下 双 圆柱 绕 流 得 到 其 升 / 阻 
丰 又 数 曲线 , 取 其 具有 代表 性 的 间距 比 2.4.10 时 的 
如 图 7~ 图 9 所 示 。 可 以 看 出 ,上 游 阻力 值 普遍 
励 间 下 游 图 柱 。 主 要 原因 是 上 游 圆 柱 对 下 游 圆 柱 起 
遮挡 的 作用 ,而 下 游 圆柱 的 升力 系数 振动 幅 值 
上 游 圆柱 的 振动 幅 值 ,上 游 圆柱 后 方 的 涡 结构 
m uu e 随 
D 的 不 断 增 加 ,上 游 圆 柱 的 涡 脱 落 逐 渐 充 分 发 
闫 落 的 尾 涡 碰 撞 到 下 游 圆 柱 从 而 对 下 游 圆柱 产 
击 作用 ,导致 了 下 游 圆柱 的 升 /阻力 值 也 和 逐渐 增 
WOM L/D 4 时 ,上 游 圆柱 升 /阻力 系数 发 生变 
化 ,曲线 变化 逐渐 平稳 。 这 主要 是 由 于 随 着 L/D 的 
逐渐 增加 ,下 游 圆柱 对 上 游 圆 柱 的 干扰 逐渐 减 小 。 
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图 7(a) 上 游 圆柱 


ChinaXiv 合 


作 期 刊 


升 /阻力 系数 


3 
时 间 /s 


图 7(b) 下 游 圆 柱 
图 7 ZD=2 时 圆柱 的 升 / 阻 力 系数 


Fig.7 Lift and drag coefficients of upstream and 


downstream cylinders with L/D - 2 
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(b) 下 游 圆 柱 
图 8 LD=4 时 上 下 游 圆 柱 的 升 / 阻 力 系数 
Fig.8 Lift and drag coefficients of upstream and 


downstream cylinders with L/D =4 
通过 数值 模拟 得 到 串联 双 圆 柱 不 同 间距 下 的 上 
下 游 时 均 阻 力 系数 值 如 图 10 所 示 。 可 以 看 出 : 当 
2«L/D«3.5 时 , 随 着 LD 的 不 断 增加 ,上 游 圆 柱 表 
现 出 渐渐 减 小 的 规律 ,而 下 游 圆柱 表现 出 先 增 天、 后 


减 小 的 规律 ; 当 ZD =4 时 ,上 游 圆 柱 和 下 游 圆柱 的 
时 均 阻 力 系 数 同时 发 生 了 较 大 的 变化 ; 当 4<LID= 
10 时 ,上 游 圆 柱 时 均 阻 力 系 数 不 再 随 着 L/D 的 增 大 
而 发 生 较 大 的 变化 ;下 游 圆 柱 在 L/D =4 时 的 时 均 
阻力 系数 为 0.31, L/D =6 时 的 时 均 阻 力 系数 为 
0. 37 ,相差 仅 为 0.06 ,总 体 未 发 生 太 大 的 波动 。 
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O 图 9 L/D =10 时 上 下 游 圆柱 的 升 /阻力 系数 


Fig.9 Lift and drag coefficients of upstream and 
downstream cylinders with L/D = 10 

图 11 ~ 图 13 分 别 为 三 维 串 联 双 圆 柱 绕 流 在 代 
表 性 的 间距 比 24、10 时 ,流向 截面 和 展 向 截面 瞬时 
三 维 涡 结构 分 布 图 。 可 以 看 出 : 当 L/D =2 时 ,由 于 
两 圆柱 间距 过 小 ,上 游 圆 柱 后 还 未 完全 形成 完整 的 
尾 涡 便 附着 在 下 游 圆 柱 上 ; 当 L/D =4 时 , 双 圆 柱 间 
距 增 加 后 上 游 圆 柱 后 方 的 尾 涡 逐 渐 充 分 发 展 并 产生 
明显 的 涡 脱 落 现 象 ,由 于 后 方圆 柱 阻 碍 了 该 尾 涡 的 
继续 发 展 ,上 游 脱 落 的 尾 涡 在 碰撞 到 下 游 圆 柱 时 会 
产生 许多 复杂 细小 涡 洲 ; 当 LD =10 时 ,在 间距 比较 
大 时 上 下 游 圆 柱 后 都 各 自发 展 出 完整 的 尾 涡 结构 , 间 


距 比 的 增 大 让 圆柱 间 的 尾 涡 结构 变 得 越 来 越 复杂 。 
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图 10 串联 双 圆 柱 时 均 阻 力 系数 


Fig.10 Time average drag coefficient of tandem double cylinders 


(b) JE fs] T 
图 11 LD=2 时 的 双 圆 柱 尾 涡 结构 分 布 
Fig.11 Wake structure distribution of 
double cylinder with L/D 22 


MT 


Ca) 流向 截面 


(b) 展 向 截面 
图 12 LD=4 时 的 双 圆 柱 尾 涡 结构 分 布 
Fig.12 Wake structure distribution of 
double cylinder with L/D =4 
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comes 


Ca) 流向 截面 


(b) 展 向 截面 
图 13 L/D =10 时 的 双 圆柱 尾 涡 结 构 分 布 
Fig.13 Wake structure distribution of 
double cylinder with L/D = 10 

738 A R SR p p — BER FD RT HE. FR E 

时 均 阻 力 系数 值 ,并 将 其 与 单 圆柱 以 及 串联 双 
圆 桂 进行 对 比 , 结 果 如 图 14 所 示 。 可 以 看 出 :串联 
绕 流 在 两 种 不 同 的 布置 方式 下 ,上 游 圆柱 的 时 
网 村 力 系数 变化 规律 较为 接近 ,说 明 在 上 游 圆 柱 后 
充 增 加 圆柱 对 其 阻力 系数 影响 不 大 。 中 游 圆 柱 在 间 
BÉEE 2 & L/D <4 时 的 时 均 阻 力 系数 呈 减 小 趋势 ,在 
AKUD <10 时 的 时 均 阻 力 系数 变化 较 小 。 这 说 明 
在 2<LD «4 时 增加 圆柱 会 使 得 中 游 圆 柱 时 均 阳 力 
ZO iit 4 e L/D «10 时 增加 圆柱 对 中 游 阻力 
系数 影响 不 大 。 下 游 圆柱 在 L/D «4 时 ,两 种 布置 
疗 下 的 时 均 阻 力 系数 变化 稍 有 差别 ;但 在 L/D >4 
时 a 两 种 布置 方式 下 的 时 均 阻 力 系数 的 变化 规律 基 
本 接近 。 这 说 明 ; 当 L/D «4 时 ,在 下 游 圆 柱 前 方 增 
加 圆柱 对 下 游 的 时 均 阻 力 系数 有 影响 ; 当 L/D >4 
时 ,在 下 游 圆 柱 前 方 增加 圆柱 对 其 时 均 阻 力 值 几乎 
没有 影响 。 
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图 14 串联 三 圆柱 时 均 阻 力 系数 


Fig.14 Time average resistance coefficient 


of three cylinders in series 
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图 15 ~ 图 17 给 出 了 三 维 串 联 三 圆柱 绕 流 在 间 
距 比 为 2.4、10 时 的 流向 截面 和 展 向 截面 瞬时 三 维 
涡 结 构 分 布 图 。 可 以 看 出 : 当 L/D =2 时 ,由 于 增加 
了 圆柱 ,间距 比较 小 时 上 游 圆柱 后 面 的 尾 涡 基本 都 
被 中 游 圆 柱 所 阻碍 ,这 使 得 中 游 圆 柱 后 方 只 有 少量 
的 尾 涡 出 现 。 当 L/D =4 时 , 随 着 间距 比 的 增加 ,中 
游 加 柱 对 上 游 圆 柱 的 尾 涡 阻碍 作用 减 小 ,中 游 圆 柱 
后 方 开 始 出 现 尾 涡 结构 。 当 L/D =10 时 ,在 较 大 的 
间距 比 下 各 个 圆柱 对 彼此 尾 涡 的 干扰 基本 可 以 忽 
略 ;各 圆柱 后 方 都 可 以 观察 到 完整 的 尾 涡 脱 落 ,每 个 
圆柱 后 方 的 尾 流 洲 涡 尺度 虽然 很 小 但 是 却 更 加 亲 
乱 ,尤其 下 游 圆 柱 后 方 尾 涡 结 构 最 为 明显 。 
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Ca) 流向 截面 


TA 


(b) JE [8] HIR 


图 15 LD=2 时 的 三 圆柱 尾 涡 结构 分 布 
Fig.15 Wake structure distribution of three cylinders with L/D =2 
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(a) 流 向 截面 


(b) 展 向 截面 
图 16 LD=4 时 的 三 圆柱 尾 涡 结构 分 布 
Fig.16 Wake structure distribution of three cylinders with L/D =4 


Ca) 流向 截面 


(b) 展 向 截面 
T= 图 17 LD=10 时 的 三 圆柱 尾 涡 结构 分 布 
> Fig.17 Wake structure distribution of three 
"T cylinders with L/D 210 
d 
45i ide 
le 


CO 通过 采用 LBM-LES 方法 对 三 维 圆柱 绕 流 进行 
模 轴 ,得 到 单 圆柱 绕 流 和 电 联 多 圆柱 绕 流 的 三 维 尾 
油 结 构 形态 及 升 /阻力 系数 曲线 ,并 进行 分 析 后 得 出 
RE, 

ND 三 维 贺 柱 绕 流 的 尾 涡 结构 比 二 维 更 直观 ,有 
SEM 


感 。 高 雷诺 数 下 单 圆 柱 绕 流 的 尾 迹 发 展会 变 得 


定 ;串联 多 圆柱 绕 流 尾 迹 存在 两 种 不 同 的 形态 ， 
HIMANEN L/D =4。 当 ZAD<4 时 ,各 
圆柱 后 方 不 会 发 生 涡 脱 落 ; 当 L/D=4 时 ,上 游 圆 柱 
百 田 开始 产生 明显 的 涡 脱 落 现象 , 且 脱 落 的 尾 涡 碰 
撞 到 后 方圆 柱 会 有 小 尺度 的 涡 产 生 。 

2) 前 方圆 柱 脱落 的 尾 涡 撞 到 后 方圆 柱 时 ,会 使 
后 方圆 柱 的 升力 系数 振幅 值 高 于 前 方圆 柱 ; 在 临界 
HEELE L/D =4 时 ,串联 多 圆柱 的 上 下 游 圆柱 时 均 
阻力 系数 会 发 生变 化 , 且 上 游 圆 柱 的 升力 系数 曲线 
也 会 变 规律 。 

3) 当 串 联 多 圆柱 绕 流 中 各 圆柱 间距 小 于 临界 间 
EEEE L/D =4 时 ,在 下 游 圆 柱 前 增加 圆柱 会 使 中 游 
圆柱 时 均 阻 力 系 数 逐 渐 减 小 ,并 且 对 下 游 圆 柱 时 均 
阻力 系数 有 影响 ;而 当 各 圆柱 间距 大 于 临界 间距 比 
L/D =4 时 ,在 下 游 圆柱 前 增加 圆柱 对 其 时 均 阻 力 系 
数值 基本 没有 影响 。 
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